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^ ■ Abstract. We propose here a model for solidification of a liquid contents of an elastic 

O ' bottle in a freezer. The main goal is to explain the occurrence of high stresses inside 

^ I the bottle. As a by-product, we derive a formula for the undercooling coefficient in 

p ^ terms of the elasticity constants, latent heat, and the phase expansion coefficient. We 

<^ ■ investigate the well-posedness of the three-dimensional model: we prove the existence 

QQ ■ and uniqueness of a solution for the corresponding initial-boundary value problem 

CN . which couples a PDE with an integrodifferential equation and an ordinary differential 
inclusion ruling the evolution of the phase parameter. Finally, we prove some results 

Oh ' on the long time behavior of solutions. 
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o : 

cn ■ 1 Introduction 

^ ■ We derive a simple model for solid-liquid phase transition of a medium inside an elastic 

0> . container. The main goal is to give a qualitative and quantitative description of the 

. . ! interaction between volume, pressure, phase, and temperature changes in the situation 

I that the specific volume of the solid phase exceeds the specific volume of the liquid 

^ ' phase. We compute the undercooling coefficient for the special case of water and ice. 

H : 

. , There is an abundant classical literature on the study of phase transition processes, 

see e.g. the monographs [3], [1], [20] and the references therein. In [5], the authors 
proposed to interpret a phase transition process in terms of a balance equation for 
macroscopic motions, and to include the possibility of voids. Well-posedness of an 
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initial-boundary value problem associated with the resulting PDE system is proved 
there. 

The microscopic approach has been pursued in [H] in the case of two different 
densities Qi and §2 for the two substances undergoing phase transitions. The evolution 
of a liquid substance, e.g., water, in a rigid container subject to freezing is described 
by a mass balance in Eulerian coordinates, an entropy balance, and a phase field 
equation. The flow is governed by a counterpart of the Darcy law. Since the density 
g2 of ice is lower than the density gi of water, experiments - for instance the freezing 
of a glass bottle filled with water - show that the water pressure increases up to the 
rupture of the bottle. When the container is not impermeable, freezing may produce a 
non-homogeneous material, for instance water ice or sorbet. This particular example 
is treated in [6] where the model is presented and a suitable variational formulation of 
the resulting nonlinear and singular PDE system is solved. In the present paper, we 
have also other applications in mind. 

Let us also mention the papers [17] and [18] dealing with macroscopic stresses 
in phase transitions models, where the different properties of the viscous (liquid) and 
elastic (solid) phases are taken into account and the coexisting viscous and elastic 
properties of the system are given a distinguished role, under the working assumption 
that they indeed influence the phase transition process. The model there includes 
inertia, viscous, and shear viscosity effects (depending on the phases), while thermal 
and phase expansion of the substance are neglected. This is reflected in the analytical 
expressions of the associated PDEs for the strain u and the phase parameter x'- the 
X-dependence, e.g., in the stress-strain relation leads to the possible degeneracy of the 
elliptic operator therein. In [17] and [H] , respectively, local existence (in the 3D case) 
and well-posedness (in the ID case) for the corresponding initial-boundary value prob- 
lems are proved. Finally, we can quote in this framework the model analyzed in [13] 
and [H], which pertains to nonlinear thermoviscoplasticity: in the one- dimensional 
(in space) case, the authors prove the global well-posedness of a PDE system, incor- 
porating both hysteresis effects and modeling phase change, which however does not 
display a degenerating character. 

Here, in Section [2], we derive a completely different model without referring to any 
microscopic balance laws, and deal exclusively with physically measurable quantities. 
We assume that the displacements are small. This enables us to state the system in 
Lagrangian coordinates. The main difference with respect to the Eulerian framework 
in [6] is that in Lagrangian coordinates, the mass conservation law is equivalent to 
the same constant mass density in liquid and in solid, but the specific volumes of 
the liquid and solid phases are different. For simplicity, we assume that the speed of 
sound, specific heat, heat conductivity, viscosity, and thermal expansion coefficient do 
not depend on the phase, the evolution is slow, and the shear viscosity, shear stresses, 
and inertia effects are negligible. The process is driven by energy balance, quasistatic 
momentum balance, and a phase dynamics equation. Still in Section [2], we verify the 
thermodynamic consistency of the model, and in Section[3]we study the equilibria. We 
observe there that a pure solid state can only be reached if the external temperature 
is below a certain threshold, which is lower than the freezing point and depends in 
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particular on the elasticity of the boundary. For water and ice, we explicitly compute 
the undercooling rate, which turns out to be around 5% if the container is rigid. For 
intermediate temperatures between freezing point and undercooling limit, there exists 
a continuum of distinct equilibria with mixtures of solid and liquid. If in this situation 
the bottle breaks, an instantaneous solidification takes place. 

The well-posedness of the three-dimensional model is investigated in Section HI 
and the asymptotic stabilization of the process is proved in Section [5l 

2 The model 

As reference state, we consider a liquid substance contained in a bounded connected 
bottle C with boundary of class C^'^ . The state variables are the absolute 
temperature 6 > 0, the displacement u e M^, and the phase variable x ^ [0? 1] • The 
value X = means solid, x = 1 means liquid, x ^ (0, 1) is a mixture of the two. 

We make the following modeling hypotheses. 

(Al) The displacements are small. Therefore, we state the problem in Lagrangian 
coordinates, in which the mass conservation is equivalent to the condition of a 
constant mass density qo> 0. 

(A2) The substance is compressible, and the speed of sound does not depend on the 
phase. 

(A3) The evolution is slow, and we neglect shear viscosity and inertia effects. 
(A4) We neglect shear stresses and gravity effects. 

In agreement with (Al), we define the strain e as an element of the space T^^^ 
of symmetric tensors by the formula 

£ = V.u:= ^(Vu+(Vu)^). (2.1) 

Let 5 G T^sym denote the Kronecker tensor. By (A4), the elasticity matrix A has the 
form 

Ae = X{e:5)S, (2.2) 

where " : " is the canonical scalar product in T^y^ , and A > is the Lame constant 
(or bulk elasticity modulus), which we assume to be independent of x by virtue of 
(A2). Note that A is related to the speed of sound vq by the formula vq = \/X/ Qq. 

We want to model the situation where the specific volume VgoUd of the solid phase 
is larger than the specific volume Vuquid of the liquid phase. Considering the liquid 
phase as the reference state, we introduce the dimensionless phase expansion coefficient 
« = {VsoUd - Viiquid)/Viiquid > 0, aud we define the phase expansion strain i by 

six) = ^{l - X)S . (2.3) 



3 



We fix positive constants Cq (specific heat), Lq (latent heat), 9^ (freezing point at stan- 
dard atmospheric pressure), 70 (phase relaxation coefficient), P (thermal expansion 
coefficient), and consider the specific free energy / in the form 

/ = coe(l-\og(^))+^{{e-~e{x)):Sr-^{e-e^)e:S (2.4) 

where / is the indicator function of the interval [0, 1] . 

To derive the balance equations, we first proceed formally, assuming that the 
temperature is positive. This assumption will be justified in the subsequent sections. 
The stress tensor a is decomposed into the sum cr" + cr^ of the viscous component 
a'" and elastic component <t^. The state functions cr^' ,cr'^ , s (specific entropy), and 



e (specific internal energy) are given by the formulas 

a-" = i^iet:S)S (2.5) 

ct' = go^ = iHe:S-a{l-x))-Pie-e,))S, (2.6) 



e = f + 9s = co9 + —{e:S-a{l-x)? + -9ce:S + Lo{x + Hx))A2-8) 

2Qo Qo 

where u > is the volume viscosity coefficient. The scalar quantity 

p := -ust :d-X{s:d-ail- x)) + P{0 " ^c) (2.9) 

is the pressure and the stress has the form cr = —pS. The process is governed by the 
balance equations 

diver = (mechanical equilibrium) (2.10) 

^oCf + divq = (T : St (energy balance) (2-11) 

~loXt ^ d^f (phase relaxation law) (2-12) 

where is the partial subdifferential with respect to x > and q is the heat ffux vector 
that we assume in the form 

q = -kV9 (2.13) 

with a constant heat conductivity n > 0. The equilibrium equation (12.101) can be 
rewritten in the form Vp = , hence 

p{x,t) 

Pstand + P{t) , (2.14) 

where Pstand is the constant standard pressure, and P is a function of time only, 
which is to be determined. We assume the external pressure in the form P^xt = 
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Pstand + Po with a constaiit deviation . The normal force acting on the boundary is 
— (cr + Pe^((5)n = (P(t) — Pq)6iy = {P(t) — pQ)n, where n denotes the unit outward 
normal vector (this notation is slightly ambiguous: the first two terms in this vector 
identity involve left multiplication of a vector by a matrix, while the last term is a 
vector multiplied by a scalar). We assume an elastic response of the boundary, and a 
heat transfer proportional to the inner and outer temperature difference. On dfl , we 
thus prescribe boundary conditions for u and 6 in the form 

(P(t)-po)n = k(x)u, (2.15) 
q-n = h{x){e-er) (2.16) 

with a given symmetric positive definite matrix k (elasticity of the boundary), a 
positive function h (heat transfer coefficient), and a constant > (external tem- 
perature). This enables us to find an explicit relation between divu and P. Indeed, 
on dQ we have by (12.151) that u ■ n = (P(t) — po)]<.~^ {x)n[x) ■ n{x) . Assuming that 
k^^n ■ n belongs to L^{dQ) , we set 

^ k"i(x)n(x) ■n(x)ds(x), (2.17) 



Kr Jan 
and obtain by Gauss' Theorem that 

Unit):= [ divu(x,t)dx=^(P(t)-po). (2.18) 

Under the small strain hypothesis, the function divu describes the local relative vol- 
ume increment. Hence, Eq. (12.181) establishes a linear relation between the total rela- 
tive volume increment t/n(t) and the relative pressure P{t)—pQ. We have £ : 6 = divu, 
and thus the mechanical equilibrium equation (12.141) . due to (12.91) and fl2.18p . reads 

z/divut + A(divu - a{l - x)) - Pid - 9^) = -po - KrUn{t) . (2.19) 

As a consequence of (12.41) . the energy balance and the phase relaxation equation in 
(l2ll|l -( l2A2l l have the form 

QqCqOi — i%/S.6 = u^divuf)'^ — pOdivUf — [a\{divu — a{l — x)) + QoLo)xt ^ (2.20) 
-QoloXt e aX{divu-a{l-x)) + QoLo(l--^ + dI{x)] , (2-21) 



where d denotes the sub differential. For simplicity, we now set 

c:=goCo, l-=golo, L:=QqLq. (2.22) 

The system now completely decouples. For the unknown functions 9^Xi ^i-nd U = 
div u , we have a closed system of one PDE and two "ODEs" (note that mathematically, 
dl{x) is the same as Ldl{x)) 

cdt-KM = uU^ - peUt- {aX{U -a{l-x)) + L)xt, (2.23) 

i,Ut + XU = aX{l - x) + I3{e - 6,) - Po - KrUnit) , (2.24) 

-^Xt e aX{U-a{l-x)) + L(l-pl+dIix), (2.25) 



5 



with Un{t) = fg^U{x,t) ds{x) , and with boundary condition fl2.16p . fl2.13p . To find 
u, we first define $ as a solution of the Poisson equation A$ = U with the Neumann 
boundary condition V$ ■ n = KrUQ{t)]<.~^ {x)n{x) ■ n{x) . With this $, we find u as 
a solution to the problem 

divu = in nx (0, oo) , (2.26) 
u • n = , (G + V<l> - KrUnk-^n) x n = on dQ x (0, oo) , (2.27) 

and set u = u + V$. Then u satisfies a.e. in Q the equation divu = U , together 
with the boundary condition (12.151) . that is, u = i^r%k~^n on dQ. 

For the solution to fl2.26p - fl2.27p . we refer to [HI Lemma 2.2] which states that for 
each g G H^^'^{dil)^ satisfying Jg^g - nds{x) = there exists a function u G H^{^1)^ , 
unique up to an additive function v from the set V of divergence- free H^{Q) functions 
vanishing on dQ , such that div u = 0inf2, u = gon dQ . In terms of the system 
fl^r^ - fl^^TT]) . it suffices to set g = ((V$ - A'rf/f^k^^n) x n) x n and use the identity 
(b X n) X n = (b ■ n) n — b for every vector b . Moreover, the estimate 

mf \\u + ^r\\H^n) < C \\g\\m/2(9n) < C\mHHn) (2.28) 

holds with some constants C, C . The required regularity is available here by virtue 
of the assumption that is of class C^'^ , provided belongs to H^^'^{dVL). Note 
that a weaker formulation of problem (12.260 - 02.270 can be found in [H Section 4]. 

Due to our hypotheses (A3), (A4), we thus lose any control on possible volume 
preserving turbulences v G V . This, however, has no infiuence on the system (12.230 - 
fl2.25p . which is the subject of our interest here. Inequality (12.281) shows that \ eV 
can be chosen in such a way that hypothesis (Al) is not violated. 

In terms of the new variables 6,U,x, the energy e and entropy s can be written as 
e = c^9 + ^{U-a{l-x)f + -9cU + Lo{x + Hx)), (2-29) 

s = co\og(l)+^x+-U. (2.30) 
The energy functional has to be supplemented with the boundary energy term 

Erit) = + (2-31) 

The energy and entropy balance equations now read 

f goe{x,t)dx + Er{t)\ = [ h{x){9r - 9) ds{x) , (2.32) 



/ Jan 

Qost + div^ = -L^ + lx'^ + -Uf > 0, (2.33) 
A f g^s{x,t)dx = I ^{9r-d)ds{x) (2.34) 

ft;|V6'P 7 2 ^ rr2\ 1 
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The entropy balance (I2.33P says that the entropy production on the right hand side is 
nonnegative in agreement with the second principle of thermodynamics. The system 
is not closed, and the energy supply through the boundary is given by the right hand 
side of ( Km . 

We prescribe the initial conditions 

e{x,o) = e\x) (2.35) 

f/(x,0) = f/°(x) (2.36) 
X{x,0) = x\x) (2.37) 

for X & Q, and compute from f l2.29l) -( l2.30l) the corresponding initial values e° , i?p , and 
s° for specific energy, boundary energy, and entropy, respectively. Let = QqC^ dx , 
= QqS^ dx denote the total initial energy and entropy, respectively. From the 
energy end entropy balance equations (12.321) . (12.341) . we derive the following crucial 
(formal for the moment) balance equation for the "extended" energy gQ{e — 6rs) : 

J (eg + ^{U - a{l - x)? + .BOM + Lx) {x, *) dx + ^ (ua(t) + (2,38) 

+ / / ^i9r-9Yix,r)dsix)dT 
Jo Jdn 9 

= E' + E^-erS' + er J^(^c\og(^pj+jX + PUj {x,t)dx. 

We have log(^/^e) = log(^/2^r) -log(^c/2^r) < (^/2^r) - 1 -log(^e/2^r) , hence there 
exists a constant C > independent of t such that for all t > we have 

Jje + u') (x, t) dx + J^J^ + f + ^) ^2.39) 
+ f [ ^ier-e)\x,T)ds{x)dT < c. 

Jo Jdn ^ 



3 Equilibria 

It follows from (I2.16P and (I2.23P that the only possible equilibrium temperature is 
9 = 9t , and the equihbrium configurations Uoo, Xoo for U, x satisfy for a.e. x E fl the 
equations 

Af/oo(x)-aA(l-Xoo(x)) = /3{er-e,)-po-Kr [ U^{x')dx' , (3.1) 

Jn 

-Af/oo(x) + aA(l-Xoo(x)) G -(l-^]+dI{xoo{x)), (3.2) 
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as a consequence of fl2.24p . 02.251) . hence 

^(j-l^-PiOr-0,)+po + Krj^Uoo{x')(^x' e dI{xoc{x)) a.e. (3.3) 
The equihbrium pressure Poo is given by (12.181) . that is, 

Poc=Po + Kr I f/oo(x') dx' . (3.4) 
Jn 

Integrating Eq. (13. ip over Vt yields 

{\ + KrM) [ Uooix') dx' = \n\{/3ier - Oc) - po) + a\ [ (1 - Xoo(a;')) dx' . (3.5) 
Jn Jn 

Hence, a necessary and sufficient condition for Xoo{x) to be an equihbrium phase 
distribution reads 

(3.6) 

Let us introduce a positive dimensionless parameter 

_ a'XKr\n\ 

- L{X + Kr\Q\) - ^ ' 

Assume ffist that /3/(A + A'r|f2|) and Po/{\ + AV|f2|) are neghgible with respect to 
the other terms. We then rewrite Eq. (13.61) in a simpler form 

1^-1+ |^^(l-Xoo(x'))dx' G dI{xoo{x)) a.e. (3.8) 

We distinguish three cases: 



er>9c 



Then (13. 8p can only be satisfied if Xoo = 1 a-e., hence, by (13.50 . Uoo = a.e., 
and by (13. 4p . the pressure Poo is in equilibrium with the external pressure. We 
only have the liquid phase in Q and the system is stress-free. 



< 1 and < (1 - d)e^ 



Then, similarly, (13. 8p can only be satisfied if Xoo = a.e., hence 

aX\n\ 



Uc^(x') dx' = , Uoo{x) = a- 



oo 

n 





n 


aX 


X + K 


r 


n 


X + Kr 





We only have the solid phase subject to a balance between a positive volume 
expansion Uoo and pressure Poo — Po = -^r|^|f^oo- 

In the limit case Kr (stress-free boundary condition, i.e. infinitely soft 
bottle), we get Poo Po, Uoo — > a, d 0. Hence, a measures indeed the 
relative volume expansion in the stress-free case. Similarly, in the limit case 
Kr — > oo (rigid bottle), we have Poo — po ^ cxX, Uoo — > 0. In this case, aX is 
the pressure difference between inside and outside the bottle. 
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(1 - d)9, <9r<9c 

Set (i* = l — {9r/0c) < d. Then every function Xoo with values in [0, 1] satisfying 
the condition j^{l — Xoo{,x')) dx' = d^^/d is an equihbrium. Hence, in this 

temperature range, we have a large number of possible equilibria. 

We thus observe stable undercooled mushy regions in a nonzero temperature 
range, and full solidification only takes place if the temperature is below the value 
(1 — d)6c. Theoretically, we cannot exclude the case d > 1, which would mean that 
the solid phase can never be achieved. We show now that in the case of water and 
ice, which is relevant for applications, the undercooling coefficient d is less than 1. 
Approximate values of the physical constants are listed in Table [Tj see [7] . 

The maximum of d is achieved in a rigid bottle (i.e. oo). By Table [T] 

we have a = {V,ce - Vwater) jV^ater = 0.09, A ^ 2.25 ■ 10^ J/m^. Using Eq. (12:221) we 
obtain L = qqLo ^ 0.33 ■ 10^ J/m^, hence d = a^\/L ^ 5.5%. Note that the standard 
atmospheric pressure is about 10^ J/m^, while the pressure inside the bottle attains 
aA ~ 2 ■ 10* J /rrt" . This corresponds to a mass of 20 kilograms pressing by gravity on 
each square millimeter. 



Specific volume of water 


Vwater — ^ / Qo 


10-3 


m? /kg 
m? /kg 


Specific volume of ice 




1.09 ■ 10-3 


Speed of sound 


vo = \/\/ go 


1.5 ■ 103 


m/ s 


Freezing point 




273 


K 


Specific heat 


Co 


4.2 ■ 103 


J/{kgK) 


Latent heat 


Lo 


3.3-105 


J/kg 


Thermal expansion coefficient 


p/x 


2.0 ■ 10-^ 





Table 1: Physical constants for water 



In reality, some values of the constants are different in water and in ice (the 
specific heat, for instance, is only 2 ■ 10^ J/ {kg K) in the ice). A phase field model 
without mechanical effects for this situation was considered in [12]. Also the speed 
of sound in ice is about the double of the one in water. We can in principle state 
the problem with coefficients depending on 6 and x here, too, but this would lead 
to serious technical difficulties that we want to avoid here. Moreover, in water and 
ice, the thermal expansion coefficient (3 is not constant and depends strongly on the 
temperature as well as on the phase. It may even become negative for temperatures in 
a right neighborhood of the freezing point. The values given in Table [T] are obtained 
by a rough linearization in order to have an idea about the orders of magnitude. 

For the coefficient P we compute the estimate P/qo = Vq[3/\ ~ 450 J/ {kg K) , 
while Lo/{a6c) = L/{QQa9c) ~ I'i'iOO J / {kg K) . Let us define a new constant : = 
Lo — j3a6c/ Qo as a small (3.3%) correction to the latent heat Lq. Using (13. 4p . we may 
rewrite (13.31) as 

^ (^^ - 1^ + Poo e dl{xoo{^)) a.e. (3.9) 
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We now show that (13. 9p contains the Clausius-Clapeyron equation, cf. P, Book 5, 
Chapter 5] or [Jlj pp. 124-126]. The pressure P^o is defined as the difference 5P 
between the absolute pressure and the standard pressure. The phase transition takes 
place at temperature 9^ if the right hand side of (13. 9p vanishes. The temperature 
difference is 66 = 6r — Oc, and we get the Clausius-Clapeyron relation in the form of 
Eq. (288) of [SP, that is, 

SP QoL/s Lg 



59 aOc Oc{V^ater " Vic 



(3.10) 



For general (3 > and po, we have an analogous classification as above. We 
introduce further dimensionless quantities 

^~ L{\ + Kr\n\)' ~ L{X + Kr\n\)' ^ ^ 

The counterpart of (13.81) reads 

(1 - /3) (^^ - 1^ + c.^ + p ^(1 - Xooix')) dx' e dl{xoo{x)) . (3.12) 

Assuming that /3 < 1, we thus observe pure liquid for 9r > Oc{l — uj/{l —/?)), while 
pure solid corresponds to 6^ < 6c{l — {u + d)/{l — (3)). The dimensionless external 
pressure deviation uj can be assumed small. However, /3 is a material constant, and 
the condition (3 <1 might be restrictive. Again, for water and ice, the maximal value 
(3 = {Lq — Lp)/Lq ^ 0.033 corresponding to Ky = shows that the infiuence of 
thermal expansion on the undercooling coefficient is negligible. 

4 Existence and uniqueness of solutions 

We construct the solution of (I2.24p - (l2.25p by the Banach contraction argument. The 
method of proof is independent of the actual values of the material constants, and we 
choose for simplicity 

L = 2, c = 9c = a = P = -f = K = \ = u = l. (4.1) 

System (I2.23I) - (I2.25I) with boundary condition (I2.16P then reads 

J 9tw{x)dx + j Ve-Vw{x)dx = J (^U^ - OUt - (^U + x + l)xt) w{x) dx {4:.2) 

h{x){e - er)w{x) ds{x), 



'an 

Ut + U + x + KrUn{t) = 6 - po , (4.3) 

Xt + U + x + dI{x) 3 26-1, (4.4) 

where (14.21) is to be satisfied for all test functions w G W^''^{Q) and a.e. t > 0, while 
(I4.3l) - (l4.4p are supposed to hold a.e. in Qao ■= ^ x (0, oo) . 

In this section we prove the following existence and uniqueness result. 



10 



Theorem 4.1 Let < 6^, < < 6* and po ^'M. be given constants, and let the data 
satisfy the conditions 

f/°,X° e L°^in), < x\x) < 1 a.e. 

Then there exists a unique solution {9, U, x) to (^T^-(^T^, li2.35\) - [2.3'7^ , such that 
9 > a.e., X e [0,1] a.e., U,Ut,Xt,9,l/9 E L°°in^), 9t, A9 E L'^{Vt^), and V9 E 
L^{0,T;L\n))nL\n^). 



Remark 4.2 For existence and uniqueness alone, we might aUow the external tem- 
perature 9r to depend on x and t, and assume only that it belongs to the space 
^loc (0' L'^i.dVt)) n L'^^iydVt X (0, oo)). For the global bounds, the assumption that 
be constant plays a substantial role. 



The proof of Theorem 14. II will be carried out in the following subsections. Notice 
first that the term U'l — 9Ut — (f^ + X + ^)Xt on the right hand side of fl4.2p can be 
rewritten alternatively, using (14.41) and (14. 3p . as 

U^-9Ut-{U + x + l)xt = U^-9Ut + Xt-'20xt (4.5) 

= -{x + U + po + KrUn)Ut-{U + x + l)Xt, 
We now fix some constant R > and construct the solution for the truncated system 

/ 9twix)dx+ [ V^-Vw(x)dx = / {Uf + Xt -QRid)iUt + 2xt))w{x)dx (A.G) 
Jn Jn Jn 

-I h{x){9 - 9r)w{x) ds{x) Vu; G iy^'^(l]) , 
Jan 

Ut + U + x + KrUnit) = Qr{9)-Po, (4.7) 

Xt + U + x + dI{x) ^ '2Qr{9)-1 (4.8) 

first in a bounded domain Qt := x (0,T) for any given T > 0, where Qr is the 
cutoff function Qr{z) = minj^"*", i?} . We then derive upper and lower bounds for 
9 independent of R and T , so that the local solution of (14. 6 1) - (14. 8 p is also a global 
solution of (1T^ - (10|) if R is sufficiently large. 



4.1 A gradient flow 

In a separable Hilbert space H with norm | ■ | , consider a gradient fiow 

v{t)+dij{v{t)) 3 fit), v{0) = v\ (4.9) 

where : H ^ [0,oo] is a. proper convex lower semicontinuous functional such that 
lim|^l^oo "^("y) = +00, dip is its sub differential, and E Dom-?/;, / E L^(0, oo;i7) 
are given. A classical existence and uniqueness result in [21 Theoreme 3.6] states that 
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for every T > there exists a unique solution v G C{[0,T]; H) to (14.91) such that 
V E L2(0,r; H), and 

We prove here the following Lemma. 

Lemma 4.3 Let f,f belong to L^{0, oo; H) . Then limt^oov{t) = 0. 



Proof. For each h > and a.e. t > we have 



1 d 

2 dt 



v{t + h)~ v{t) 



hence 



h 



< 



\V[S] 



fit + h)- fit) 



v{t + h) - v{t) 



h 



< 2 



|/(r)|Kr)|dr 



(4.10) 



for almost all < s < t. Hence, the function t ^ 2 /q |/(T)||t'(T)| dr — |?)(t)P is 
almost everywhere equal to a nondecreasing function in (0,cx3). We are thus in the 
situation of [T5l Proposition 5.2], which gives the desired statement. ■ 



V 



We apply the above result to the case H = L^(f2) x L^(f2) , and 

U 
X 

in v2 



(f/ + X - 1)' + (f/ + 2x)(l - ^r) + Hx) ) dx 



(4.11) 
(4.12) 



Kr 



f 



Qr{0) - 9r 
2{QR{e) - 9r] 



(4.13) 



where is a suitable constant such that ip{v) > for all v , and 6' is a given function. 
The initial condition v'^ is given by (I2.36p . (12.371) . We have 



c 



G dipiv) 



7] = U + x-Or + Kr J^Udx + po, 
C G U + x + l-2er + dI{x), 



(4.14) 



and we see that Eqs. (14. 71) - (14. 8 p with 9 replaced by 9 can be equivalently written as 
a gradient flow (14.91) . (I4.1ip -( l4.13p . For its solutions, we prove the following result. 



Proposition 4.4 Let the hypotheses of Theorem \4.1\ hold, and let a function 9 G 
Lf^^{0,oo;L'^{Q)) be given. Let {U,x) be the solution of ^.11\ )- ^T^ . Then 

there exists a constant Cq, independent of x,t and R, such that a.e. in fioo we have 



\Uix,t)\ + \Utix,t)\ + \xt{x,t)\ < Coil + R) . 



(4.15) 



12 



Let furthermore 61,62 G L^q^(0, 00; he two functions, and let (f/i, Xi); (^2, X2) 

he the corresponding solutions of ^4-9^ , ^.ll^ - l^.l^ . Then the differences 6^ = 
61 — 62, Ud = Ui — U2 , Xd = Xi ~ X2 satisfy for every t > and a.e. x E Vt the 
inequality 

j^md)t\ + \{Xd)t\){x,T)dT<Co{l+t) (\Ux,T)\+t\6d{T)\2) dr, (4.16) 
where the symhol | ■ I2 stands for the norm in L'^{fl) . 



In what follows, we denote by Ci, C2, ■ ■ ■ any constant independent of x, t and R. 

Proof. Put Xfi^t) = — x{x' ,t)) dx' . Integrating (14.71) with 6 replaced by 6 over 
fl yields 



Un + {1 + Kr\n\)Un = Xn + (Q/?(^) - 1) dx - poM a.e. 

Jn 

Since x attains values in [0, 1] , we easily obtain + < Ci{l + R) a.e. Equation 
(14. 7p now has a right hand side bounded by a multiple of 1 + -R, hence \Ut \ + \U\ < 
6*2(1 + R) a.e. To obtain the same bound for \xt\ , it suffices to multiply ( 14. 8 p by Xt- 
This completes the proof of (I4.15p . 

To prove ( I4.16p . we rewrite ( 14. 9p . ( I4.1ip -( l4.13p as two scalar gradient flows 

Ut + dMU) = a, (4.17) 
Xt + dMx) ^ (4.18) 

where MU) = \U\ ^2 = \x^+I{x), a = QR{6)-x-Po-KrUn, b = 2Qr{6)-1-U. 
Consider now two different inputs. As above, we denote the differences {}i — {}2 by 
{}d for all symbols {}. By [ini Theorem 1.12], we have for all t > and a.e. x G f2 
that ^ ^ 

[ md)t\ + \{Xd)t\){x,T)dr<2 [ (|arf| + |6,|)(x,r)dr. (4.19) 
Jo Jo 

We multiply the difference of (I4.17P by Ud, the difference of (I4.18P by Xd, and sum 
them up to obtain that 

iUd)tUd + iXd)tXd + iUd + Xd)^ + KrUndUd< \6d\i\Ud\ + 2\xd\) a.e. (4.20) 

We first integrate (I4.20p over Q. Using the symbol | ■ I2 for the norm in L'^{Q), we 
get for a.e. t > that 

{\Ud\l + \Xd\l) + KrUL < \6d\2i\Ud\2 + 2\xd\2) < V5\6d\2 {\Ud\l + \Xd\lf' ■ 

(4.21) 

Hence, -;^(|t/d|2 + IXdli)"*^^^ — V^l^dh a.e., and integrating over t, we find that 

{\Ud\l + \Xd?^"\t)<V^ f\6d{r)\2dr. (4.22) 

Jo 
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This implies in particular that 

\Und{t)\ < VWA [ \0d{r)\2dr. (4.23) 
Using again fl4.2Up . we find for a.e. {x,t) G Qoo the inequality 

2^ {\Ud\' + IXdl') {x,t) < (^Kr\Und{t)\ + {\Ud\ + 2\xd\){x,t) . (4.24) 

This is for almost all x G an inequality of the form (d/ dt)(Y'^{t)) < 2c(t)Y{t), 
y(0) = 0, with Y = (|?7rf|2 + which implies Y{t) < c(r) dr for alU > 0. 

Hence, 

(|f/dr + |Xd|')'^'(a:,t) <Ci£(|^,(x,r)|+t|^^rf(r)|2) dr a.e. (4.25) 
This enables us to estimate the right hand side of fl4.19p and obtain the bound 

j\\(Ud)t\ + \{xM{x,T)dT < C^j^ ((1 +t)|^rf(x,r)| +t(l + t)|^,(r)|2) dr (4.26) 
for G and all t > 0. This completes the proof. ■ 

4.2 Existence of solutions for the truncated problem 

We construct the solution of f l4.6p -( l4l8l) for every i? > by the Banach contraction 
argument on a fixed time interval (0, T) . 

Lemma 4.5 Let the hypotheses of Theorem \4.1\ hold, and let T > and R > 
be given. Then there exists a unique solution {6, U, x) to ^.6\ )- ^T^ , l{2.35\) - [2.37^ , 
such that U G W^'°°{nT) , 9 > a.e., Xt, 0,1/9 G i:°°(fir), Ot, A9 G ^^(fir), and 
Ve G L°°(0,T; L'^{n)). 

Proof. Let 9 G L'^{Qt) be a given function, and consider the system 

[ 9tw{x)dx+ [ V9 ■Vw{x)dx = [ {U^ + Xt -QBi9){Ut + 2xt))w{x)dx (4.27) 
Jo, Jn Jn 

-[ h{x){9-9r)w{x)ds{x) G Vr^'^(fi) , 
Jan 

Ut + U + x + KrUn{t) = Qr{9)-po, (4.28) 
Xt + U + x + dI{x) ^ 2Qr{9)-1. (4.29) 

Equations fl4.28p -( l4.29p are solved as a gradient fiow problem from Subsection 14.11 
while fl4.27p is a simple linear parabolic equation for 9 . Testing fl4.27p by 9t , we 
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obtain by Proposition 14.41 that 

[ [ e^dxdt + supess ( [ \V9\^dx + [ h{x){9 - 9rY ds{x)] (4.30) 
Jo Jn te{o,T) \Jn Jan J 

< T\n\ (Co(l + i?)(2Co(l + R) + 3R)f =: Mr . 

Hence, we can define the mapping that with 6 associates the solution 6 of fl4.27|] - 
fl4.29p with initial conditions f l2.35l) -( l2.37l) . We now show that it is a contraction on 
the set 

Et,r := {§ E L^iVtr) : conditions ^i^-^M) hold} , (4.31) 

where 



Ve G L°^{0,T;L\n)); 



(4.32) 
(4.33) 



rT 




Jn 



te{o,T) \Jn 



^jMxdt + supess / \\/9\^dx+ / h{x){e - Orf ds{x) ] < Mr; (4.34) 



an 



e{x,0) = e%x) a.e. 



(4.35) 



Let 61,62 be two functions in Ht,_r, and let {6i,Ui,xi) , (6*2, f^2, X2) , be the corre- 
sponding solutions to fl4.27p - fl4.29p with the same initial conditions 6^,U^,x^- We 
see from fl4.30p that 61,62 belong to Hj'/j. Integrating Eq. fl4.27p for 61 and 62 with 
respect to time and testing their difference hj w = 6(i := 61 — 62, v/e obtain, using 
Proposition 14.41 that 

6lix,t) dx + ^ 

< C,{l + R) [ ( [\\{Ud)t\ + \{Xd)t\ + m{x,T)dT)6,{x,t)dx a.e. (4.36) 









\ 


/ 6d{x, t) dr 


dx + h{x) 


/ 6d{x, r) dr 


ds{x) j 


/o 


Jan 


Jo 





n \Jo 



From fl4.16p and Minkowski's inequality, it follows that 



(|(t/.)*| + |(x.)t|)(-,r)dr 



< 0,(1 + tf / |^",(r)|2dr 



1/2 



By Young's inequality, we rewrite (14.360 as 



6l{x,t) dx + — 



V / 6d{x,T)dT 



dx + h{x) 



an 



6d{x, r) dr 



ds(x) 



< C^{l + R^){l + tf / \dd{T)\ldT a.e. 



(4.37) 
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Set e^{t) = /o*|^d(r)|^dr, Q^t) = /q |^d(r)|i dr . Integrating fOTI) with respect to 



time, we obtain 



e2(t) < C5(l + i?') / (l + r)^e2(r)dr. (4.38) 
Jo 

We set Cr := (C5(l + R^)/Q) and introduce in L°°(0,T) the norm 

llu'lic := sup e-^«(^+")VMI- 

Then ||0||^ < |||0||^, and hence the mapping 6 9 is a contraction in L'^{Qj') with 
respect to the norm induced by || ■ \\c- The set 'Et,r is a closed subset of L'^{Qt) ■ This 
imphes the existence of a fixed point 6 G '^t,r, which is indeed a solution to (14.61) - 
(14.81) . The positive upper and lower bounds for 6 follow from the maximum principle. 
Indeed, the right hand side (14.51) of (14. 6 p is bounded from above by Cq{1 + Ry and 
from below by —^{9~^y. Let us define the functions 

e\t) = 9* + Ceil + R)H , 9\t) 



2 + 9J 

For every nonnegative test function w and a.e. t G (0,T) we have 

9tw{x)dx+ / V9 ■Vw{x)dx+ / h{x){9 - 9r)w{x) ds{x) (4.39) 



an 



< Ceil + Rf I w{x)dx, 
Jn 



9tw{x)dx+ / \/9 ■Vw{x)dx+ / h{x){9 - 9r)w{x) ds{x) (4.40) 
n Jn Jan 

> -\ I {9^fw{x)dx, 



2 



! 9\w{x)dx+ I V9^ ■Vw{x)dx+ I h{x){9^ - 9t)w{x) ds{x) (4.41) 
Jn Jn Jan 

> Ceil + RY I w{x)dx, 
Jn 

9\w{x)dx+ [ \/9^ ■Vw{x)dx + I h{x){9^ - 9t)w{x) ds{x) (4.42) 
n Jn Jan 

< -\ I {9''fw{x)dx. 



2 



n 



We now subtract (I4.4ip from (14.390 and test hy w = {9 — 9^)^ , which yields the 
pointwise bound 9{x,t) < ^"(t). Similarly, we subtract (I4.40p from (14.420 and test by 
w = {9^ — 9)'^ . We thus have the inequalities 

9\t) < 9{x,t) < 9^{t) a.e., (4.43) 

which complete the proof of Lemma 14.51 ■ 
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4.3 Proof of Theorem 14.11 



The unique solution (9, U, x) to (^-i^, fl03D - (l07l) exists globally in the whole 
domain Qoo- We now derive uniform bounds independent of t and R. Take first for 
instance any R > 26* . By fl4.43p . we know that the solution component 6 of fl4.6p -( H^ 
remains smaller than i? in a nondegenerate interval (0,T) with T > 6* /{Cq{1 + R)"^) . 
Let (0,To) be the maximal interval in which 6 is bounded by R. Then, in (0,To), 
the solution given by Lemma [4.51 is also a solution of the original problem fl4.2l) - fl4.4l) . 
Moreover, due to estimate fl2.39l) . we know that 6 admits a bound in L°°{0, Tq; L^{fl)) 
independent of -R. In order to prove that Tq = +oo if R is sufficiently large, we need 
the following variant of the Moser iteration lemma. 

Proposition 4.6 Let Q C be a bounded domain with Lipschitzian boundary. 
Given nonnegative functions h G L^{dQ) and r G L°°(0, oo; with a fixed 

q > N/2, \r\L°°{o,oo;Li{n)) ='■ r* , an initial condition G L°°{Q) , and a boundary 
datum vr G L°°{dQ x (0, oo)), consider the problem 

Vt-Av + v = r{x,t)H[v] a.e. m x (0, oo) , (4.44) 



under the assumption that there exist positive constants m, Ho,Cf,V,Vr, Eq such that 
the following holds: 

(i) The mapping Ti : L'^^iVL x (0,oo)) — > L'^^i^ ^ (0,oo)) satisfies for every v G 
L'^^{VL X (0, cxd)) and a.e. {x,t) G x (0, cxd) the inequality 



v{x,0) 



Vf ■ n 



V 



,0 



h{x) {f{x,t,v{x,t))—vr{x,t)) a.e. on 9f2 x (0, oo) , (4.45) 
' a.e. inn, (4.46) 



v{x,t)n[v]{x,t) < Ho\v{x,t)\ (l + \v{x,t) \ + [ ^{t-T)\v{x,T)\dT 



) 



where ^ G l/r^'"'^(0, oo) is a given nonnegative function such that 



m<-mm a.e. 



(4.47) 



(ii) f is a Caratheodory function on Q x (0, oo) x M such that f{x,t,v)v > Cjv"^ 
a.e. for all v 



(iii) |f°(x)| <V a.e. in Q . 



(iv) \vr{x,t)\ < Vr a.e. on dQ x (0, oo) . 




(0, oo; W''\n)) n L^^{n x (O, oo)) 
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Then there exists a positive constant C* depending only on |/i|Li{an); Cf, Hq such 
that 

K'(t)|L-(n) < C*max{l, V,\/r,£'o} for a.e. t > 0. (4.48) 

Remark 4.7 As a consequence of (lOT]) . we have ^(t) < ^(0)e-«W* for all t > 0, 
hence ^{t) dt < 1. As a typical function satisfying (14.471) . let us mention for 
example 

n 

e(t) = ^^?^5^r,e-™'=* (4.49) 

with any < mi < ■ ■ ■ < m„ and > , k = 1, . . . ,n. 

We split the proof of Proposition 14.61 into several steps. 

Lemma 4.8 Let ^ be as in Proposition \4 ■ 6[ let x E W^^^{0,oo) and y G Ll^^{0,oo) 
be nonnegative functions, and let a > , C > 0, 5 G (0, 1) be given constants. Set 
H = min{a, ^(0)(1 — S)} , and assume that for a.e. t>0 we have 

x{t) + ax{t) + y{t) < C + 5 f ^(t - r) i/(r) dr . (4.50) 

Jo 

Then x{t) < max{x(0), C//i} for all t> 0. 

Proof Set z{t) = /J^(t - r)2/(r) dr. Then (l/^(0))i(t) + z{t) < y{t) a.e., hence 

x{t) + ax{t) + ■^^^i't) + (1 - < C a.e. 

With /i as above, we have 

which yields 

<t) + -^^t) < max |x(0) + ^^(0), ^| , 
and the desired inequality follows easily. ■ 

Lemma 4.9 Let Ti be as in Proposition "JTB, and let \ ■ \p denote the norm in Lp{Q) 



for 1 < p < oo . Let v G L^^{Qoo) md p, s > 1 be arbitrary. For {x,t) G ^loo set 
hp = 7^[f] V \v\P~^ . Then, for all t > we have 



\hp{t)U < H,{^^ + ^\v{t)\l+^-j\{t 
where ^ is as in (I4.47p . 



r)\v{r)\l,dr , , 
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Proof. We have for a.e. (x, t) G Qoo that 

\hp{x, t)\<Ho {\v{x, t) + 2\v{x, t) |P + t)) 

where 

1 ^ 



Wp{x,t) 



P \Jo 



^{t-T)\v{x,T)\dT 



1/p' / rt 



l/p\P 



< - [\{t~T)\v{x,T)\^dT. 

P Jo 

Here, we have used Holder's inequahty with conjugate exponents p,p' and Remark ITTl 
The assertion now follows from Minkowski's inequality 



^{t-TM;T)\PdT 



< / at-r)Hr)\l^dT. 
Jo 



Lemma 4.10 Let the hypotheses of Proposition hold, and let \ ■ |oo,p denote the 
norm m L°°(0, oo; LP(fi)) . Let v G W^^^{0,oo; {W^'^y{n)) f] Ll^{0,oo;W^^^{n)) n 
L'^^{fl X (0, oo)) be a solution to /li4.44\ )^ [4-4(^ snc/i that \v\oo,p < oo for some p > I . 
Then there exists a constant C > independent of v and p such that 



oo,2p < {Cp'+Y'" max {1, V, Vr/Cf, \vUJ 



where 



N{q + 1) 



2q-N 



(4.51) 



(4.52) 



Proof. We test (14.441) by ^ to obtain, using Lemma [4.91 that 

— 4- f \v(x,t)\^Pdx + [ \\/\v\PWx,t)dx+ [ \v(x,t)\^Pdx (4.53) 

2p dt p^ 



Jdn 

We estimate the boundary integral using Young's inequality 

v^v\v\^p-^ < ^P:ilcf\v\^p + —c\''pv'r 

' ' - 2p ^' ' 2p f ^ 
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Set Vp{x,t) = \v{x,t)\P. Then 

-^^\vMl + 2|V^,(t)|^ + 2p\v,it)\l (4.54) 

By the Gaghardo-Nirenberg inequahty, [12], there exists a constant G such that for 
all 5 > and t > we have 

\vMl' < G{5\Vv,{t)\l + S~'\vMl) < GiSlVvMl + S-'lvlXp), (4-55) 
with b given by (14.521) . We now choose 6 such that 

Qpr*HoG6 = 1 , 

and obtain 

f^K{t)\l + \'^Mt)\l + MvMl < Cril + iVr/Gfr^+p'-^'lvlX,) (4.56) 



1 /"* 

— / at-r)\Vvp{T)\ldT a.e., 

Qp Jo 



with a constant C7 depending only on .^(0), Gf, |/i|Li(aQ), r* , Hq, and G. We 
now use Lemma ITS) with x{t) = \vp(t)\2, y{t) = |Vwp(t)|2, a = 2p, 6 = l/{6p), 
/X > min{2, (5/6)^(0)}, and G = Gr{l + {Vr/Gf^P + p^+'>\v\^J^^p) , which yields that 

\vMl < Cs{l + V'P+{Vr/GfrP + p'+'\v\Xp) (4.57) 
with a constant Cg independent of p and t, and (14.511) immediately follows. ■ 
We are now ready to finish the proof of Proposition 14.61 

Proof of Proposition \4-6\ For A; = 0, 1, 2, . . . set yk = max {l, V, Vr/Gj, |f |oo,2fc} • We 
have yo < max {1, V, Vr/Gf, Eq} , and, as a consequence of Lemma 14.101 

This yields 

log|/fc+i < logyk + 2-('=+i)(logC' + kil + h) log 2) . 

Hence, 

n 

logy„<logyo + 5^(2"'=(logC+(A;-l)(l + 6)log2)) . (4.58) 

fc=i 

The sum on the right hand side of (I4.58p is convergent, and we easily complete the 
proof. ■ 
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We now finish the proof of Theorem 14.11 by showing that Tq introduced at the 
beginning of this subsection is +00 if R is sufficiently large. In (14. 3p . set again 
U^{t) = J^U{x\t)dx'. Then 

Unit) + {I + KrM)Un{t) = [ (6 - x){x' ,t) dx' - \n\po a.e. 

Jn 

By ( 12.391) . the right hand side of this ODE is uniformly bounded independently of R, 
hence |t/n(i)| < Cg in (0,To). Using (14.31) once again, we obtain that 

\U{x,t)\ < Cio(^l + e^-'e{x,T)dT^ a.e., (4.59) 

\Ut{x,t)\ < Cii ^1 + ^(x,t) + ^ e"~*^(x,r)dr^ a.e., (4.60) 

hence also (cf. (lO]) ) 

\Xt{x,t)\ < Ci2(l + 0{x,t) + j^e^-'e{x,T)dT^ a.e. (4.61) 

As in (14.51) . we rewrite the right hand side of Eq. (14.21) as 

-{X + U + PO + KrUnWt -{U + x + l)Xt • 

By ( 12.391) . the function U is in L°°{0, oo; L'^{Q)) and the bound does not depend on R. 
Eq. ( 14. 2p . with 6 added to both the left and the right hand side, thus satisfies the 
hypotheses of Proposition 14.61 for = 3 and q = 2. This enables us to conclude that 
9{x,t) is uniformly bounded from above by a constant, independently of R, so that 
6 never reaches the value i? if i? is sufficiently large, which we wanted to prove. By 
(I4.59p - (l4.6ip . also U , Ut, and Xt are uniformly bounded by a constant. 

We proceed similarly to prove a uniform positive lower bound for 6. Set Rq := 
sup^, and in Eq. (14. 6 p with R > Rq put w = —w/9, w G W^''^{Q). For a new 
(nonnegative) variable v{x,t) := logi?o — log6'(a;, t) we obtain the equation 



Vtw{x) dx + j Vf ■ \/w{x) dx + I h{x) ( ~ 1 ) ui{x) ds{x) (4.62) 
n Jn Jan \ " 



n 



U^ + Xl \'^^\\ut + 2xt]w{x)dx. 



9^ 



We now set 



(tt2 I 2 \\7f)\^ \ 



and check that the hypotheses of Proposition are satisfied with f{v) = {6t/Ro){g" — 
1), vr = {Ro — Or)/Ro, r = 1, and vH[v] < 2C\v\, where C is a common upper bound 
for Ut and Xt- Hence, v is bounded above by some v* , which entails 6 > Roe'"" . 
This concludes the proof of Theorem 14. 1[ 
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5 Long time behavior 



In order to emphasize the relation to Section [31 we keep the original physical constants 
as in fl2.23|] - fl2.25l) . We prove the following statement. 

Proposition 5.1 Let the hypotheses of Theorem \4.1\ hold, and let the constants /3,uj 
introduced in Ii3.11\) satisfy the condition 1 — P > max{0,ci;}. Then we have 

[ ( [ [el + f/2 + xl + |V^|2) dx+ f h{x){9 - 9rf ds{x)] dt < oo ,(5.1) 
Jo \Jn Jdn J 



lim ( [ (Uf + xt + \^d\'^){x,t)dx+ [ h{x){e-erfix,t)ds{x)] 
\Jn Jdn J 



0. (5.2) 



Furthermore, the functions Un{t) = J^U{x,t)dx, X^{t) = J^(l — x(x, t)) dx converge 
to their equilibrium values as t ^ oo . 

The function on the left hand side of (15.21) is almost everywhere equal to a function 
of bounded variation. The limit is to be understood in this sense. 

We see in particular that the temperature converges strongly in W^''^{fl) to its 
equilibrium value 6r , and the total ice contents Xq as well as the pressure converge 
to a constant as t ^ oo. For temperatures above the freezing point or below the 
undercooling limit, this means, in agreement with the discussion in Sectional that also 
both xi^jt) ^^d U{x,t) converge strongly in L^{Q) (hence, strongly in every L^{Q) 
for j9 < oo) to their respective equilibrium values. For intermediate temperatures, 
only the limit total ice contents can be identified, but we are not able to decide about 
the convergence of the individual trajectories to some of the equilibria. 

Proof. By (I2.39p . since 6 is uniformly bounded from above, we have 

/ / (f/^^ + xt + |V^|^) (x,t)dxdt + / / h{x){e -erf{x,t)ds{x)dt <oo. 
Jo Jn Jo Jdn 

(5.3) 

We rewrite Eq. fl2.23p in the form 

c6t - kM = uiUtf - peUt + JXt - ^Oxt . (5.4) 

Due to the uniform L°° upper bounds for 6 and Ut , we can test Eq. (15.41) by 9t , and 
obtain 

/ e^{x,t)dx + ^( [ \veWx,t)dx+ [ h{x){e - Orf ds{x)] (5.5) 

Jn dt \Jq Jg^ J 

< Ci3 / {U^ + Xt) dx a.e. , 
Jn 

with Ci3 independent of t, and (15. 3p with (15.51) together with (THl Lemma 3.1] yield 

lim / |V^P(x,t)dx+ / h{x){e-erf{x,t)ds{x)=0. (5.6) 
Jn Jan 
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System fl2.24l) - fl2.25l) can be again considered as a the gradient flow of tlie form fl4.9l) . 
with V and ipi^v) analogous to fl4.1ip - (l4.12p . more precisely 



uU 



(5.7) 



- ^^2 , /r,, , pn TT\ f-i 



^(^) = i^-^iU - a{l - x)y + (Lx + pe^U) (^1 - / ) + I{X) ) dx (5.8) 



^ [{L/e^)ie-er) ) 



(5.9) 



We have / G L\0, oo; L^{Q) x L'^{n)) by (El, and / G L\0,oo;L\n) x L'^{Q)) by 
(15.51) . From Lemma [4.31 we conclude that 

lim / (U^ + Xt) (x, t)dx = 0. (5.10) 

Set 6't7(t) = J^6{x,t) dx . The equation for f/^ now reads 

uUn + (A + Ar|fi|)t/o = aXX^ + - {po + Wl^^l , (5.11) 

hence 

lim ((A + KrM)Unit) - aXXn{t)) = (/3(^r - ^c) - Po)M ■ (5.12) 
From (I2.24P and 02.251) we obtain for a.e. (x, t) G fioo that 

A([/ -a{l-x)) = -yUt + /?(^ - e,) - Po - A^rf/n , (5.13) 

-7Xt e a(-i/f/i+/3(^r-^c)-Po-A'rt/f,)+L (l - ^) + - {9-9r) + dI{x) ■ 

(5.14) 

We define an auxiliary function 

A{x,t) := --fXt{x,t) + auUt{x,t)- (ap-^]{9{x,t)-er) (5.15) 



+ aKrUn{t) xTl^\ ' ^^'^^^ 

With the notation (13.71) . (13. lip we rewrite fl5.14p in the form 

^-A{x, t) + T^Xnit) + (l-p)(^^-l]+u;e dl{x{x, t)) a.e. , (5.17) 



as an evolution counterpart of the equilibrium condition (13.120 . The above computa- 
tions show that limj^oo |^(^)|2 = 0. We now prove the following imphcations: 



23 



(i) If > ^c(l -^/(l then J^\l - x{x,t)\dx as t^oo; 

(ii) If < ^c(l - iLU + d)/{l-(3)) then J^x{x,t)dx as t ^ oo; 

(iii) If e^{l -{uj + d)/{l - §)) <er < e,{l - uj/{l - /3)) then 
Xn{t) ^ m/d){{l - /3)(1 - (er/e,)) -uj) ast^oo. 

The corresponding convergence of U then follows from f l5.12p - fl5.13p . 

To prove the above statements (i)-(iii), set 

Eq. (15.171) reads 

A*{x,t) + ^^Xn{t)-x* ^dl{x{x,t)) a.e. , (5.18) 

that is, 

(^A*{x, t) + ^^Xn{t) - x*^ (X - Xix, t)) < a.e. Vx G [0, 1] . (5.19) 
Integrating over Q, we obtain for every x ^ [0, 1] and a.e. t > that 

(^l^Xn(t) - X*) (i^Xn(t)-(l-x)) <-|^^A*(x,t)(x-x(a:,t))dx. (5.20) 
The right hand side of fl5.20p tends to as t tends to oo . Hence, 

hmsup ( p^n(t) - X*) ( p^f^W - (1 - X)) < Vx G [0, 1] . (5.21) 

(i) We have x* ^ 0. The assertion follows if we put x = 1 in fl5.2ip . 

(ii) We have x* ^ 1- The argument of (i) applies if we put x = in (15.211) . 

(iii) Here, we have < x* < 1 , and it suffices to put x = 1 — x* ■ 

U 

Note that in all cases the difference U — a{l — x) converges in L^(f2) to its 
equilibrium value as t — oo. The problem if x{x,t) and U{x,t) separately converge 
in the case (iii) is still open. 
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